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SUMMARY 


This report gives a model, an iterative calculation method, and results for plane 
collisionless sheaths. These sheaths separate thermally ionized plasmas (Saha) from 
emitters of electrons (Richardson-Dushman), ions, and atoms (Saha-Langmuir) and 
change emission with their fields (Schottky). Still the model gives identical absolute po- 
tential changes through the sheath for ions and electrons, contrary to usual practice. 

For cesium, the emitter and plasma ion and atom temperatures range from 1600*^ to 
2400° K with plasma electron temperatures up to 300° K higher. Plasma electron num- 
ber densities are 10 to 10 per cubic centimeter. At these conditions work functions 

from 1. 5 to 5 volts produce absolute potential drops from 0 to 1. 5 volts for electron and 
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positive-ion sheaths. And some currents greater than 10 amperes per square centime- 
ter result. 

For this study certain generalizations apply: Sheath fields affect emissions and 
sheath characteristics significantly. The overall sheath characteristics correlate better 
with an emission Debye length (6. 9 times the square root of the emitter temperature 
(°K) divided by the plasma electron number density (cm“ )) than with the conventional 
Debye length. Where conditions comply with the assumptions of the model, effective 
sheath widths decrease from 2. 6 to 2 emission Debye lengths as overall sheath drops 
grow from less than 0. 1 to 1. 5 volts. 


INTRODUCTION 


In gaseous electronics an electrode exchanges flows of electrons, ions, and atoms 
with an ionized gas. Depending on the properties of the system, the potential of the 
solid surface generally runs above or below that of the plasma. The region between the 
plasma and the emitter builds up an excess of electrons or positive ions. Under certain 



conditions these sheaths submit to analyses based on existing theories. For example, 
relations for equilibrium ionization, emission, acceleration, and potential combine to 
form a simple picture of a collisionless adaptive region between an electrode and an ion- 
ized gas. This sheath model with its iterative calculation method and some results for 
cesium make up the present report. 

Although a sheath and a net current generally attend an emitter in gaseous electron- 
ics, they vanish at a special state of equilibrium, the Saha-Langmuir null point (ref. 1). 
Here the emitter injects electrons, ions, and atoms into the plasma at rates that exactly 
balance the departures of these particles from the equilibriated thermally ionized gas. 

In other words, the Richards on- Dushman emission of electrons (refs. 2 to 4) and the 
Saha-Lai^muir emissions of ions and atoms (ref. 1) equal the respective random parti- 
cle currents to the emitter from the Saha plasma (ref. 5). With the ionization potential, 
pressure, and matched plasma and electrode temperatures all fixed, equilibrium dic- 
tates the work function at the Saha-Langmuir null point. 

If the work function is less than the value that yields no net current or sheath, fewer 
positive ions and more electrons escape the emitter. Then an excess of negative charge 
gathers near the face of the metal and diminishes toward neutrality in the plasma, form- 
ing an electron sheath. When the work function is greater than that of the Saha- 
Langmuir nvill point, emission increases for ions and decreases for electrons. Now a 
surplus of positive charge collects at the solid surface and diminishes toward zero in the 
plasma, producing an ion sheath. In both positive-ion and electron sheaths, though, no 
net currents result if the emitter remains at the temperature of the equilibrium plasma. 

But sheaths sometimes confront electrode surfaces with great charge concentrations 
and high fields that create problems not anticipated in the Saha-Langmuir model. And if 
the system is not isothermal, if net charge flows, if conditions depart from equilibrium 
in any way, both the Saha-Langmuir and Saha theories fail - technically. Practically, 
however, these relations serve well for many nonequilibrium conditions (refs. 6 to 8). 
For certain modes of operation the Saha and Saha-Langmuir equations adequately de- 
scribe power producing thermal-plasma devices like thermionic converters (refs. 6 
to 9). 

Because electrodes as well as plasmas contribute charged particles in the therm- 
ionic diode, both of these active elements affect charge concentrations in the sheaths of 
this energy converter. Although many previous studies analyzed sheath problems (refs. 

9 to 27), those that dealt with plasmas and probes used nonemitting or nonionizing elec- 
trodes. And those that treated emitters included no ionization mechanism within the 
plasma and no increases in emission caused by the effective lowering of the work func- 
tion by electric fields (ref. 28, Schottky). 

These Schottky depressions of potential barriers at emitters often affect charge 
transport greatly in thermionic diodes. For example, a 1-volt sheath of a plasma with 
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10^^^ electrons per cubic centimeter at 3000° K probably rises about 10*^ volts per centi- 
meter near the electrode. This field cuts almost 0. 04 volt from the effective work func- 
tion and frees 25 percent more electrons than escape with no field. Such reduced bar- 
rier potentials and increased emission currents demonstrate that sheath models for 
thermionic diodes need Schottky modifications. 

This work describes a field- corrected emitter (Richardson-Dushman, Saha- 
Langmuir, Schottlq^) connected through a collisionless sheath (electron or positive-ion) 
to a thermally ionized plasma (Saha). The sheath produces identical overall changes in 
absolute potential for electrons and ions, contrary to the results of usual Schottky anal- 
yses. To illustrate the present model, plots and tables in this report give charge and 
current densities, electric fields and potentials, and widths of sheaths for cesium plas- 
mas. Electrode, ion, and atom temperatures range from 1600° to 2400° K; electron 
temperature varies from 1700° to 2700° K; work function runs from 1. 5 to 5. 0 volts: 

1 Q 1 C ' 

and electron number density goes from 10 to 10 per cubic centimeter. 

The results presented reflect the restrictions of the model: The sheath is essen- 
tially collisionless. Emitter, ion, atom, and plasma temperatures are similar. Net 
currents are small compared with random circulations. And the plasma pressures for 
cesium are below its vapor pressures at the electrode temperatures. 

These computations convert to another plasma merely by replacing the minimum 
mean free path, vapor pressure, and ionization potential of cesium with those of the 
desired element. 


THEORY 


As stated in the IN TRODUCTION , this model joins a thermally ionized plasma to a 
field- corrected emitter with a collisionless sheath. The theory begins, however, with 
the exact Saha-Langmuir solution that yields no sheath and no net current. 


SAHA-LANGMUIR NULL POINT 

The plasma model synthesized for the present analysis contains only electrons and 
singly ionized and ground-state atoms. This simple model serves well for low tempera- 
tures. In equilibrated cesium plasmas, for example, excited atoms are an appreciable 
part of the neutrals only above 3000° K. At equilibrium, the proposed composition fol- 
lows Saha’s expression for the law of mass action: 
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Because these plasma electrons (eq. (1)) form a Maxwellian distribution, one-half 
of them moves in any one direction with the average velocity (2/cTg at any 

given time. A random current density results: 
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When this flow cancels the electron emission from a surface at the plasma temperature, 
the electrode and the ionized gas equilibrate: 
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With its denominator of 2, the Saha-Langmuir relation (eq. (3)) applies when the 
statistic weight of the ion is one -half that of the atom. This situation occurs in ground- 
state, singly ionized alkali-metal vapors. Here the emissions described by equation (3) 
exactly replace ion and atom losses from the ionized gas. Because no sheath intervenes, 
a given plasma demands one particular work function for true equilibrium with an elec- 
trode. Work functions that produce these Saha-Langmuir null points for some cesium 
plasmas appear in figure 1. 
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SHEATHS 


Generally, the work function of the emitter fails to satisfy the requirements for a 
Saha-Langmuir null point. Then a sheath adapts the plasma to the electrode and equal- 
izes the countercurrents in the isothermal system. Sheaths usually occur in hetero- 
thermal combinations of emitters and ionized gases also. But, in the strictest sense, 
the present development applies only where all temperatures are alike. 

For these calculations charge balance and zero electric field prevail in the bulk 
ionized gas. Zero reasonably approximates the fields needed for the comparatively 
small current densities in the plasma. As the distance from the electrode increases, 
the sheath conditions approach those of the plasma. In practice, though, various par- 
ticle and group interactions preclude this asymptote, and much of the sheath potential 
drop occupies a relatively small distance. This rationale leads to the concept of an ef- 
fective sheath width with an approximate charge balance at the plasma boundary. Fur- 
ther discussion later will clarify this idea. 


Virtual Schottky Emitter 

The present study analyzes a sheath attached to an electrode that emits electrons, 
ions, and atoms. The sheath lowers the emission barrier for the charge it accelerates 
and frees more accelerated particles than escape from a sheathless emitter. A 
positive-ion sheath releases more electrons; an electron sheath looses more positive 
ions. Consideration of this effect led to the use of a new electrode model called a vir- 
tual Schottky emitter. On its surface, particle potentials differ in absolute value from 
those in the plasma by the sheath drop for the accelerated charge. But this is a decep- 
tively simple definition. Perhaps the virtual Schottky emitter requires additional ex- 
planation, beginning with a description of traditional theoretic thermionic emission. 

Without external fields, an electron converts kinetic energy to potential while leav- 
ing a conductor. The total energy change for an exodic electron equals the work func- 
tion of the emitter. Schottky suggested that the electron works against the attraction of 
its mirror image to escape. The electron and its image are equally but oppositely dis- 
tant from the conductor surface. As a result, the charge pair attract each other with a 
force inversely proportional to four times the square of the distance of the electron from 
the emitter, e /4x'^. Of course this force has no meanii^ at the surface. Furthermore, 
the image effect applies only for a single charged particle departing from the plane face 
of a conductor, precluding other particles and external fields. But experiments verify 
image forces in the presence of other particles and of external fields (ref. 6). 

Thus, when an electron moves toward a conductor from very far away, it loses po- 
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tential energy equal to the image force times its distance from the surface, e /4x. Be- 
cause this change lowers the free-space potential energy of the electron, it subtracts 
from the work function, e<p - e /4x. Dividii^ this potential energy by the electronic 
charge yields the motive, which figure 2(a) shows. The mirror image effect requires 
the use of the motive concept because in this situation the usual definition of potential 
has no meaning. 

K a charged particle moves in superposed electrostatic fields, its charge divided 
into the potential energy resulting from the mirror image and all external fields is its 
motive. But here, as throughout thermionics, ’’potential” is also used to mean 
’’motive.” The negative gradient of any motive equals the electrostatic force on the 
particle divided by its charge. A particle field of this type results from the mirror im- 
age, e/4x^. 

When an external field acts on an emitter, it adds vectorially to the image force. 
With a constant applied field that accelerates electrons, the emission motives become 
cp - (e/4x) - Ex for the electrons and (p - (e/4x) + Ex for the ions; these functions ap- 
pear in fi^re 2(b). On the lower curve the maximum <p - (eE)^/^ occurs at a distance 
e/2(eE)^' ^ from the conductor surface. This depression of the emission barrier, ef- 
fectively reducing the work function for the electron, is the Schottky effect. 

If the applied field accelerates the positive ion, the upper curve of figure 2(b) de- 
scribes the electronic potential. Then, as the lower curve indicates, the emission bar- 
rier for positive ions decreases. 

Traditionally, thermionics treats the work function and the Schottky effect as sur- 
face, not space, phenomena. This assumption leads to motives that terminate at the 
emitter face like those in figure 2(c). In the present model the electron and ion poten- 
tials in the plasma lie at equal displacements above and below the origin. For usual 
analyses, though, emission barriers are the work function for the ion and the field- 
corrected work function for the electron. When the plasma potentials subtract from 

those at the emitter, the ion- sheath drops become cp^ - cp for ions and 
1/2 ^o ^ 

<p - (p^- (eEg) ' for electrons. Thus, particles with equal but opposite charges un- 
dergo different absolute potential changes while passing through the same sheath. Fur- 
thermore, the traditional emitter of figure 2(c) produces net currents through the sheath 
in an isothermal system. Obviously then, the usual thermionic model fails in a sheath 
analysis. 

A motive like that of figure 2(b) seems necessary to combine the image and sheath 
potentials, as the diagrams of figure 2(d) show. For these motives no net charge flows 
when all temperatures are alike. But now the emission mechanism fuses complexly and 
inseparably with the sheath transport mechanism, preventii^ good independent estimates 
of emitter characteristics. The transition region looks like a single sheath for one 
charge and a double sheath for the other, with the image potential extending through the 
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sheath and the plasma. This means that plasma calculations include effects of surface 
potentials for all electrodes and containers. 

The extension of the image potential through the sheath and the plasma, however, 
lacks both theoretic and experimental support. As reference 9 puts it, ”. . . the 
mirror-image force is of little or no consequence at a distance from the conducting sur- 

O 

face of 10“ meter. . . ." Because the superposition of figure 2(d) is unreal and com- 

plicated, it seems a poor replacement for the traditional model of figure 2(c), which is 
incorrect but simple. 

Thermionics needs an effective emitter that separates emission from sheath and 
plasma processes yet submits to the second law of thermodynamics. The virtual 
Schottky emitter offers the compromise shown in figure 2(e). In this model, electrons 
and positive ions pass through the same absolute potential change across the sheath, 

AVq - (eEg)^/^ for the ion sheath and AV^ - (-eEg)^/^ for the electron sheath. Si- 
multaneous emissions of electrons, ions, and atoms encounter the same effective work 

I/O 1/2 

functions, <p - (eEg) ' for the ion sheath and <p + (-eEg) ' for the electron sheath. 
The sheath potentials for ions and electrons at the emitter occur at slightly different 
distances from the solid surface. But all of the other emission models for electrodes in 
plasmas also have this problem. Finally, no net charge flows for isothermal systems 
involving virtual Schottky emitter. 

The virtual Schottky emitter, a theoretic construct that works for sheath calcula- 
tions, is the basis of the present analysis. 


Parameters 

Calculations for both the positive-ion and electron sheaths begin with the specifica- 
tion of certain variables: emitter work function and temperature Tg; plasma electron, 
ion, and atom temperatures T^ p, Tj p, and p; plasma electron (and ion) number 
density p = Nj p; and the ionization potential I, minimum mean free path A, and 
vapor pressure at the electrode temperature (for a maximum arrival rate) p^p for the 
plasma chemical. Values of these variables selected for the present solutions for 
cesium were 


= 1. 5 to 5 V 

Te = Ti p = Ta p = 1600°, 2000°, 2400° K 
Tg p(^Tj p) = 1700° to 2700° K 
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in nil nil I 


N. p = 10^^, 10^^, 10^^ 10^® cm"^ 
e, ir 


1 = 3. 898 V 


X = 0 (10%“^^pjcm 


Pvp = 


(ref. 29) 

anti log^Q(+10. 71914 - 0. 51978 log^gT^, - 3920. 38/Tg) 


133. 322 


torr (ref. 30) 


Current Densities 


A complete set of parameters leads ultimately to the determination of sheath fields 
and potentials from charge densities contributed by the currents. Random circulations 
within the plasmas come directly from the selected values of the variables. Because 
sheath conditions influence the virtual Schottky emitter, though, current densities of 
emissions require iterative calculations. The equations used are presented in the fol- 
lowing sections. 

Positive-ion sheath. - The plasma random current densities for the electrons, ions, 
and atoms (equivalent random current density from the Saha eq. (1)) are, respectively. 
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( 6 ) 


The following electron, ion, and atom current densities for the virtual Schottky 
emitter derive from the field- corrected Richardson- Dushman and Saha-Langmuir equa- 
tions. Although the expressions for electron emission need modification for fields only, 
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those for ion and atom emission require adjustments in the cesium arrival rates for the 
sheath drops. In addition to the flow of neutrals from the plasma, the total cesium sup- 
ply to the emitter includes ion streams that the sheaths cut off or reflect. That this 
adaptation further complicates the Saha- Langmuir relation seems obvious from a com- 
parison of equation (3) with equations (8), (9), (16), (17), (19), (20), (24), and (25). The 
current density equations are as follows: 
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(9) 


Ion sheaths yield the following electron, ion, and overall net current densities; 
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J = ie + ii 


(12) 


When all the temperatures are alike, the net current densities reduce to zero for 
electrons, ions, and atoms: 
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Then expressions for isothermal ionization ratios result: 
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Electron sheath. - Equations (4) to (6) give the random plasma current densities for 
the electron sheath as well as for the ion sheath. 

But the current densities from the virtual Schottky emitter for electrons, ions, and 
atoms change: 
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( 20 ) 
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The electron, ion, and overall net current densities are 
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(23) 


When all the temperatures are alike, equations (13), (14), and (15) apply for the 
electron, ion, and atom, respectively. New expressions for isothermal ionization 
ratios result, however: 
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Charge Densities 


With the preceding current densities for the plasma and the virtual Schottky emitter, 
adding transmitted and reflected flows divided by appropriate average velocities yields 
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charge densities. The average velocities result from integrations of (1) a half- 
Maxwellian velocity distribution for the retarded particles, (2) a haK-MaxweUian trun- 
cated at the potential rise to the top of the sheath for reflected flow, and (3) a haJf- 
Maxwellian displaced by the potential drop for the accelerated charges. 

Positive-ion sheath (electron potentials are positive numbers. ) - 
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7rmj 


The two positive terms give the electron charge density at AV while the two negative 
terms represent the ions. In the plasma, AV = 0 and equation (26) becomes 
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Pp= 0 = 



( 27 ) 


The fact that the last term equals half of the ion charge density in the plasma permits a 
reasonably direct estimate of the overall sheath drop AVg„ 

Electron sheath (electron potentials are negative numbers. ) - 
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For AV = 0, 



(29) 


Here, the first term equals half of the electron charge density in the plasma and allows 
a simple estimate of AVg. 


Sheath Configurations 

Expressions for charge densities and Poisson's equation lead to sheath fields and 
potentials as functions of distance. 

Electric fields: 
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Distances: 


X = 0 at AV = AVg 


X = Xp at AV - 0 


Xp - x = 



E 


14 



‘'AVg 


The program evaluates these integrals as sums of trapezoids based on 20 equal in- 
crements of AV from zero to ^Vg. Equation (30) causes no problem, and as far as 
the machine is concerned, neither does equation (31). But when AV =0, E = 0, and 
equation (31) is undefined. The machine, however, merely takes half of the value of E 
for AVg/20. That average E used for the first AV increment produces a finite Xp. 
This distance is an effective sheath width, which includes much of the sheath structure. 
The applicability of this definition derives from the extent to which it represents a fully 
developed sheath. If the charges nearly balance at Xp, the effective sheath width offers 
utility - particularly when Xp is a few Debye lengths rather than an incessant infinity 
(refs. 31 to 34). 

Locations for these plasma, sheath interfaces result from a simple process: First, 
select the retarded flow from the emitter (ions for ion sheaths and electrons for electron 
sheaths). Next, set the number density of this flow at the effective sheath, plasma inter- 
face equal to half of the number density of the plasma electrons. Then, compute the 
sheath potential drop (eqs. (27) and (29)). Finally, use the overall voltage change in the 
field and distance integrals to determine an effective sheath thickness. 

Justification of this concept appears in DISCUSSION OF RESULTS. 


Calculation Method 

Before the main calculations a simple screening obviates many permutations of con- 
ditions precluded by the model: First, the electrode temperature must be low enough to 
maintain an effective solid emitter. Next, the cesium pressure in the plasma must be 
below its vapor pressure (ref. 30) at the emitter temperature to prevent continuous con- 
densation of cesium. And, finally, the Debye length must be shorter than the mean free 
path for cesium charge exchange (ref. 29) to make the sheath essentially collisionless. 
The last condition was selected because the distances between cesium charge exchanges 
are less than the lengths between elastic encounters of electrons and neutrals, the reso- 
nance exchanges of free and bound electrons, and the coulombic collisions at the re- 
sulting plasma conditions. The model admits no collective interactions and allows only 
collisions beginning and ending with singly ionized and ground-state atoms and electrons. 

If a set of parameters passes this screening, that system enters the primary com- 
puting process. The iterative sheath solution follows an almost classic numeric recipe: 
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(1) Solve for current densities deleting effects of sheath fields (eqs. (4) and (5) with 
(7), (8), and (27) or (18), (19), and (29)). 

(2) From the first step also obtain an estimate of the potential drop across the 
sheath (eq. (27) or (29)). 

(3) Divide the overall voltage change into 20 increments and compute charge densi- 
ties for increasing potentials through the sheath (eq. (26) or (28)). 

(4) Integrate numerically the charge density function (from the third step) over the 
voltage to yield the sheath field (eq. (30)). 

(5) Use the field at the virtual Schottky emitter and the potential fall across the 
sheath to produce new current densities for emissions (first step). 

(6) Repeat the cycle until no boundary current density that the sheath affects (eqs. 
(10), (11), (21), and (22)) changes by more than 0. 1 percent of the smallest boundary 
current density. 

(7) Integrate numerically the negative reciprocal of the electric field over voltage to 
determine sheath distances (eq. (31)). 


DISCUSSION OF RESULTS 

For all of the temperatures (1000° to 3000° K) and plasma densities (10^^ to 
10 cm ), equilibria between emitters and ionized gases occur at work functions dic- 
tated by the Saha-Langmuir theory. As figure 1 indicates, the null-point work functions 
vary directly with minus the logarithm of electron number density and almost directly 
with temperature. 

These null-point work functions also serve as plasma potentials (AV = 0) relative to 
emitter Fermi levels for all systems analyzed here except the positive-ion sheaths with 
the plasma electron temperature elevated above that of the emitter, ions, and atoms. 
For these nonequilibrium positive-ion sheaths, figure 1 shows that plasma potentials 
branch away from the Saha-Langmuir curves. The offshoots represent the variation of 
electron potentials with electron temperature when atom, ion, and emitter temperatures 
are fixed at the value for no elevation of plasma electron temperature (AT p = 0). 

Only the potential and temperature of the plasma electrons change along any one of the 
branch curves. 

For these nonequilibrium ion sheaths, potentials of the plasma electrons, with 
given values of temperature elevation and number density lie above the Saha-Lar^muir 
value at lower temperatures and below it at higher temperatures. This trend results 
from changes in ion emission, which strongly influences the formation of the positive- 
ion sheath. The changes involve two contesting effects: If the temperature for the ions, 
atoms, and emitter drops below the temperature of the electrons, random currents sup- 
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plying cesium from the plasma to the emitter decrease, and the fraction of cesium leav- 
ing the emitter as ions increases. 

At the reduced ion and atom temperature and a given charge number density, the 
plasma atom concentration is identical with that for the same system with all tempera- 
tures at the electron value. Being proportional to the square root of the temperature, 
though, the average random cesium velocity in any direction is less. Cesium ions and 
atoms move more slowly toward the emitter, lowering ion emission proportionately. 

This reduction is relatively large for a specified elevation of electron temperature at 
low temperatures. Thus, the plasma electron potential for the nonequilibrium positive- 
ion sheath rises above the Saha- Langmuir value as the general temperature level falls. 

This influence decreases as the general temperature level rises; then the second 
mechanism gains control. Decreasing the emitter temperature increases the fractional 
ion emission (eqs. (3) and (8)), often almost exponentially for the nonequilibrium positive- 
ion sheath. At high temperatures, therefore, reducing ion, atom, and emitter tempera- 
tures together yields a net gain in ion production. And the plasma electron potential 
drops below the isothermal value. 

Although elevatir^ plasma electron temperatures affects plasma potentials for the 
nonequilibrium ion sheath, it fails to change plasma potentials for the nonequilibrium 
electron sheath. This effect results from the major influence of emitted electrons on 
the structure of the electron sheath. The barrier reflects the high number densities of 
low-energy emitted electrons back to .the emitter; only an exponentially cutoff flow es- 
capes. But this small exodic current nearly balances the injection of electrons into the 
sheath from the plasma. Plasma electrons, therefore, affect these electron sheaths in- 
significantly. And the rate of escape of electrons from the emitter remains the same - 
unaltered by plasma and sheath conditions. Thus, elevation of the plasma electron tem- 
perature causes practically no change in the electron sheath at an emitter. 

For equilibria the criteria for collisionless ion and electron sheaths were satisfied 

with 10^^, 10^^, and 10^*^ electrons per cubic centimeter at 2400° K and with 10^^ and 

10^^ at 2000° K. Figure 3 describes these systems, which include work functions from 

5 4 

1. 5 to 5. 0 volts, Debye lengths from 3. 4x10 to 3. 4x10 centimeter, and plasma po- 
tentials from 3. 0 to 4. 2 volts. The parameters appear in the insert. In figure 3 the de- 
pendent and independent variables are dimensionless. The ordinates are the effective 
sheath width (xgAj)), sheath drop (AVg/AV^ = AVg/(^ - (p^)), and emitter field 
(EgXjy^(^ - <Pq)). All abscissae are Richardson-Dushman (R-D) drops 
(|e AVQ|//cT = e|<p - cp^\/kT). 

Similar results for some nonequilibrium ion and electron sheaths appear in figure 4. 
The elevations of electron temperatures account for the increased ranges of parameters 
that satisfy the collisionless criteria (1600°, 2000°, and 2400° K and 10^^, 10^^, 10^^, 
10^^ electrons/cm^). 
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The present method yields effective sheath widths for the retarded particles from 
the emitter (for electrons in an electron sheath or for ions in an ion sheath). In con- 
trast with those for equilibria (fig. 3), the effective sheath widths in figure 4 scatter 
more widely about their mean. But the trends and levels of results are comparable for 
equilibrium and nonequilibrium ion and electron sheaths. 

This paper defines a sheath system as near- equilibrium when the net current of 
each kind of particle falls below 6 percent of its corresponding random plasma current. 
Figure 5 shows the portion of the data of figure 4 that satisfy this criterion. The same 
set of curves correlates both electrons and ion sheaths, although the expanded sheath- 
width scale accentuates separations. 

The correlations of figures 4 and 5 exemplify traditional treatments of sheath char- 
acteristics - in terms of number density, temperature, and Debye lengths for plasma 
electrons. For these Saha plasmas contiguous with Saha-Langmuir emitters, however, 
emission dominates the sheaths (eqs. (13) and (14)). Where the ionized gas meets the 
electrode, the energy of the emitted particles influences charge separation to a greater 
extent than the plasma- electron temperature. The sheath data, therefore, should orga- 
nize more readily around functions of plasma- electron number density and emitter tem- 
perature. One such function is the emission Debye length « 6. 9(Tg/Ng p)^^^ 

Figure 6 shows the nonequilibrium sheath data correlated with emitter temperature 
rather than plasma-electron temperature. Results for both ion and electron sheaths 
compact around the dashed curves for equilibria above R-D drops of 2-| to 3. Below 
this region, though, reducing overall voltages moves nonequilibrium sheaths away from 
equilibria and renders these calculations inapplicable. 

The good correlations in figure 6 indicate that sheath characteristics depend 
strongly on emitter temperature, plasma-electron number density, and emission Debye 
length when emitted particles control transitions between plasmas and electrodes. 

Figure 7 shows trends of the ratio of the minority- charge density at the effective 
sheath width to that in the neutral plasma as functions of R-D drop. These calculations 
indicate that the charge imbalance decreases below 1 percent of the number density of 
plasma electrons when the R-D drop increases above 2-| to 3-2. For equilibrium and 
near -equilibrium systems, then, the effective sheath width seems a meaningful concept. 

The number density of the majority charge at the emitter appears as a function of 
R-D drop in figure 8, which includes all equilibrium and nonequilibrium data for ion and 
electron sheaths. The rise in majority charge and consequently gas density near the 
emitter results from retarded and reflected emission. Approaching the emitter through 
the sheath, the density of the minority charge falls below that in the plasma. The atom 
density remains unchanged. But the density of the dominant charge climbs to a maxi- 
mum at the emitter even though the contribution to that density by flow from the plasma 
decreases because of acceleration. 
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Figure 8 reveals that the majority -charge density at the emitter depends only on 
plasma charge number density and the R-D drop containii^ the emitter temperature. 

This simple plot further supports the concept of the emission Debye length 
(Xde ^ ^ correlation parameter. 

The trends in figure 8 indicate that the charge number density can gain two orders 
of magnitude while moving from the plasma to the emitter through a 1-volt sheath. Of 
course, at Saha- Langmuir null points, with no sheaths, the gas density remains con- 
stant throughout. 

The present results reinforce many past findings that properties based on bulk-gas 
distributions are not meaningful within a few mean free paths of a confinii^ wall 
(refs. 35 to 38). This generalization usually holds even when gradients are small at the 
boundary. Consider a simple thermally randomized gas in a container: Unless the wall 
injects the same kinds of particles in the same proportions at the same rates with the 
same energy distributions as those it receives from the gas, the distribution changes 
abruptly within several collision lengths from the wall. If the gas moves relative to the 
container, the problem is further complicated. A solid surface cannot supply random 
distributions displaced around drift velocities unless it moves at those velocities. In 
the absolute sense, then, even the mean free path, computed using the bulk-gas distribu- 
tion, means nothing at the usual confining boundary. 

The exception occurs when the boundary returns particle flows that maintain the 
distributions in the bulk-gas. In a plasma this ideal mating of a solid surface and a 
thermally randomized gas happens at the Saha- Langmuir null point. When emissions of 
electrons, ions, and atoms fail to balance their random plasma currents, however, 
sheaths generally grow. And as figure 8 implies, distributions can change rapidly even 
in static isothermal systems. 

In this study, though, the screening calculations used the distributions of the plas- 
mas to estimate collisions in the sheaths. On this basis cesium charge exchange domi- 
nates particle encounters over the ranges of plasma conditions used. To be more cor- 
rect, the coulombic tests should have utilized mean free paths averaged over the charge- 
distribution changes of the sheaths. But this test requires the prior complete solution 
of the sheath. The additional complication was deemed unwarranted for these approxi- 
mations. 

All of the previous discussion dealt with sheath stnacture. In addition to connecting 
the plasma and surface potentials, though, the sheath influences charge flow. To illus- 
trate some transport effects, figure 9 gives net current densities as functions of R-D 
drops for all nonequilibrium sheaths in this study (Tg p > Tg = T^^ P ~ ^a P^’ Because 
the elevated temperatures of plasma electrons unbalance these systems, the net elec- 
tron currents flow to the electrodes as figure 9 indicates. For this reason electron 
temperature rather than emitter temperature might serve better in the R-D drop. The 
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change tightens the correlation a bit but not enough to justify the inconsistency of repre- 
sentation. 

The curves for current densities through electron sheaths extend horizontally down 
to R-D drops near 1; below that they fall off slightly. These practically constant cur- 
rents almost equal the increases of plasma- electron random currents caused by the 
higher plasma electron temperatures. When the higher electron currents from the plas- 
mas enter the sheaths built primarily by electron emission, they accelerate toward the 
emitter, which completely absorbs them. 

Electron sheaths with small overall potential drops, where electron currents from 
the emitters and plasmas are similar, react to plasma changes. At low R-D drops the 
net negative current through the electron sheaths from the plasma falls off because of 
rise in the net ion current to the emitter. If only plasma electron temperatures in- 
crease at constant charge number densities, random currents of atoms diminish in these 
ionized gases (eq. (1)). The decreased supplies of atoms reduce ion emissions and al- 
low higher net ion currents toward the emitters. 

Figures 5 and 6 indicate that near-equilibrium exists for R-D drops greater than 

to 3. For these R-D values, the applicable current densities of figure 9 run from 
-4 2 

less than 10 to more than 10 amperes per square centimeter. 

To this point, only overall results received attention. Now figures 10 to 13 present 

detailed structures for some ion and electron sheaths separating emitters at 2000° K 

13 

from plasmas with 10 electrons per cubic centimeter. In the ionized gas, ion and atom 
temperatures are also 2000° K while electron temperatures are either 2000° or 2200° K. 
For each of the four sets of conditions, four work functions cause sheath drops to range 
from 0 to over 0. 95 volt. Tables give incremental potential changes, charge densities, 
electric fields, and distances. Tabulations of overall sheath characteristics also ap- 
pear. Finally, plots compare Richardson- Dushman and Schottky results to indicate sig- 
nificant errors caused by omitting effects of fields on emissions. Appendix B defines 
the symbols for the IBM output sheets. 

As expected, figures 10 to 13 show that absolute rates of change of electric fields 
and potentials, as functions of sheath distance, maximize at the emitter and minimize at 
the plasma. Curves for majority- charge densities bow strongly toward the origins. But 
curves for minority -charge densities bulge away from the origins for overall sheath 
drops near 0 volt and toward the origins for overall sheath drops near 1 volt. Apparently 
simple models like exponential variations of both charges with dimensionless voltages 
fail to describe this kind of sheath. 
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CONCLUDING REMARKS 


In these sheaths, which separate Saha plasmas from Schottky Saha-Langmuir emit- 
ters, emissions control charge densities. For ion sheaths raising only the electron 
temperature in the plasma with a given charge number density reduces the sheath drops. 
This reduction results from decreasing the supply of atoms to the emitter and thus di- 
minishing ion emission. For sizeable electron sheaths, however, raising the plasma 
electron temperatures has a negligible effect because the electron emission remains un- 
changed. Since the emitter dominates the sheaths, the emission Debye length 
6. 9(Tg/Ng correlates their characteristics better than the usual plasma Debye 

length6.9(T^ p/N^ p)l/2. 

Effective widths for equilibrium and near- equilibrium sheaths in the present study 
lie between 2 and 2. 6 emission Debye lengths. Corresponding values of work function 
minus plasma potential both divided by the voltage equivalent of emitter temperature 
range from near 0 to 10. As this dimensionless voltage variable closely approaches 0, 
the effective sheath width and the overall sheath voltage differences drop toward 0. For 
values of the dimensionless voltage variable greater than 2, the sheath drops stabilize 
between 0. 8 and 1. 0 times the work function minus the plasma potential. 

These generalizations apply to equilibrium and near- equilibrium results. In the 
present work, near-equilibria occur at values of the dimensionless voltage variable above 
2^ to 3. This lower limit begins the region where net particle current densities through 
the sheaths are less than 6 percent of their random circulations in the plasmas. Also for 
voltage variables above 2^ to 3, charge imbalances at effective sheath widths are below 
1 percent of plasma electron number densities. 

Accumulations of cesium in positive-ion sheaths at high work-function surfaces com- 
plicate these systems. When emitters adsorb cesium, work functions generally decrease; 
cesium arrivals from the plasma alone lower the work function of the bare emitter. But 
with positive-ion sheaths, the greater ion concentration near the electrode increases 
cesium transport and changes the work function even more. Although these phenomena 
caused no problems in the present parametric study, in practice the work function de- 
pends on the emitter, the plasma, and the sheath. 

The near- equilibrium sheaths of this study pass net electron currents to the elec- 
trodes rather than to the plasma. Some of the current densities exceed 100 amperes per 
square centimeter. At present, though, net electron flows from the emitter to the 
plasma hold more interest, particularly in thermionic energy conversion. Applying an 
appropriate field lowers the sheath potential and allows electron emission to overrun the 
excess currents caused by elevated plasma electron temperatures. Then the problems 
of accelerated charge distributions in the plasmas grow. Streams of excess electrons 
falling through sizeable potential drops enhance local ionization and often start collec- 
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tive oscillations (ref. 39), Such complications far exceed those of the present calcula- 
tions. 

But this sheath model allows net electron flows from, as well as to, the electrodes. 
At the conditions analyzed here, for example, lowered rather than elevated plasma elec- 
tron temperatures produce electronic emitters rather than collectors. 

The model fits some heterothermal systems particularly well: For retarded plasma 
particles, an increased rise in sheath potential results in smaller current exchanges with 
the electrode; then the model tolerates larger temperature differences between the op- 
posing streams. The atoms and accelerated particles from the ionized gas undergo es- 
sentially complete replacements at the plasma, sheath interface by flows from the emit- 
ter. Thus, the model requires nearly isothermal trades of these species. Such precau- 
tions aim at preserving the specified state of the plasma. Restrictions to small net 
transport densities relative to random plasma circulations assure minor perturbations. 

Within its limitations, this sheath method analyzes equilibria, estimates near- 
eqmlibria, and gives some insight into certain systems capable of high power outputs. 

But the near- equilibrium criterion limits the current densities and temperature differ- 
ences of potentially high-power systems in analyses of this kind. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 20, 1967, 

120-33-02-01-22. 
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APPENDIX A 
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T 

V 

AV 


X 

K 

X 




SYMBOLS 


electric field 

electronic unit charge 

Planck’s constant 

ionization potential 

overall net current density 

particle current density or net 
particle current density 

particle mass 

particle number density 

pressure 

absolute temperature 
potential 

potential relative to plasma 
potential 

distance from emitter 

Boltzmann constant 

mean free path for cesium charge 
exchange in emitter tempera- 
ture range 

plasma Debye length, 

»6. 9(T^ p/N^ cm 


"DE 


P 

<P 

(P 

<P 

<P 


00 


emission Debye length, 

=6. 9(Tj,/N^ p)V2 


cm 


charge density 
work function 

effective work function for virtual 
Schottky emitter 

work function (plasma potential) 
that precludes an ion sheath 

work function (plasma potential) 
that precludes an electron sheath 

Subscripts: 


a atom 

E emitter 

e electron 

i ion 

o to work function (from plasma 

potential) 

P plasma 

S overall sheath 


vp vapor pressure 
AV at AV 
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APPENDIX B 


NOMENCLATURE FOR IBM OUTPUT 


FORTRAN 

symbol 

Algebraic 

symbol 

Description 

Units 

I 

I 

ionization potential for plasma 
atoms 

V 

TE 

Te 

emitter temperature 


PHI 

<P 

work function 

V 

NEP 

P ^ ^i, P 

plasma electron number density 

cm~^ 

TEP 

^e,P 

plasma electron temperature 


TIP 

^i, P = p for these 

solutions 

plasma ion temperature 


LAMBDA 

p/^e, 

plasma Debye length 

cm 

PV 

^vp 

vapor pressure of plasma 

c 

element at Tj. (N/m^ 

torr 

Vl 33 . 322; 

LAMBDA(TE) 

^DE ~ ^(^E^^e, p)^^^ 

emission Debye length 

cm 

DV 

AV 

sheath potential measured from 
plasma electron potential 

V 

ND(DV) 

^e, A V ■ ^i, AV 

net number density of charge 
at AV 

cm"^ 

NE(DV) 

^e,AV 

electron number density 
at AV 

cm"^ 

NI(DV) 

AV 

ion number density at AV 

cm~^ 

E(DV) 

®AV 

electron electrostatic field 
at AV 

V/cm 

X(DV) 

^AV 

distance from emitter to AV 

cm 

JEE 

^e, E 

emitted electron current density 

A/cm^ 

JEP 

^e, P 

plasma electron random 
current density 

A/cm^ 

JIP 

P 

plasma ion random current 
density 

K/cn? 
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FORTRAN 

Algebraic 

symbol 

symbol 

JAP 

k, P 

J 

J 

PP 

Pp 

JIE 

^i,E 

JAE 

^a, P 

JA 

ja 

JI 


JE 

^e 

JA/JAP 

P 

JE/JEP 

^e^^e, P 

JI/JIP 

i/ p 

DVS 


XDVS 

^AVg "" ^P ^ ^ 

NAP 

''a,P 

XDLAM 

XgAn 

SC 


PHZ 

fo 

EDVS 

®AVg " ®E 

DVSRD 

AVq = <p - cp^ 

DVS/RD 

AVg/AVg = AVg/((p - 

ELM/RD 



Description Units 

2 

plasma atom equivalent random A/ cm 

current density 

2 

net current density through A/ cm 

sheath 

plasma pressure torr 

(N/mVl33. 322) 
2 

emitted ion current density A/ cm 

2 

emitted equivalent atom current A/ cm 
density 

2 

net equivalent atom current A/ cm 

density 

2 

net ion current density A/ cm 

2 

net electron current density A/ cm 


overall sheath voltage drop V 

effective sheath thickness cm 

plasma atom number density cm”' 


Schottky depression of work V 

function 

plasma potential (work function V 

for no sheath) 

electrostatic field at the V/cm 

emitter 

Richardson- Dushman overall V 

sheath voltage drop 
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FORTRAN 

symbol 

PHZZ 

DRD/KT 

NTP 

NCE 

NTE 

RD/KTE 

X/LMTE 

ELT/RD 

NEPA 


Algebraic Description Units 

symbol 

<Pqq - for equilibrium plasma potential at equilibrium V 

and electron sheaths (work fvmction for no sheath 

and no net current) 

e\cp - cp^\/kT^ 

3 

total particle number density cm 

in plasma 

-3 

total charge number density at cm 

emitter 

-3 

total particle number density at cm 

emitter 

~ I'^'^E 

Xg/Xjjg 

-3 

plasma electron number density cm 

from sheath calculations 
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Ion sheaths with plasma elec- 
tron temperatures raised 
ATg p above those of the 
ion^ and atoms represented 
by ATg p = 0 on each branch 
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Figure L - Work functions that preclude sheaths between 
emitters and cesium Saha plasmas (plasma potentials). 
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(b) Emission motives in constant external electrostatic field. 
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(c) Traditional thermionic emitter with positive-ion sheath. 
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(e) Virtual Schottky emitters. 


Figure 2. - Development of virtual Schottky emitters. 
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Figure 3 - Equilibrium (isothermal) collisionless ion (^>V>qq) and electron 
(<P<(Pqq) sheaths. 
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(a) Ion {(P<(Pq) sheaths. 

Figure 4. - Nonequilibrium (Tp p > Tr = T, p = T. p) collisionless sheaths 
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Figure 4. - Concluded. 
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Figure 6. - Nonequilibrium (Tg p > = Tj p = Tg p) collisionless sheaths (related to 

emission Debye length). 
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Figure 6. - Concluded. 
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Figure 7. - Examples of charge imbalances at nonequilibrium cesium plasma sheath 
interfaces. 
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Figure 8. - Emitter charge densities in cesium plasmas. (All 
equilibrium and nonequilibrium data for ion and electron 
sheaths are included. ) 
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Figure 10. - Plane ion sheaths between emitters and cesium plasmas at equilibrium. 
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Figure 11, - Plane ion sheath between an emitter and a cesium plasma near equilibrium. 
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Figure 12, - Plane electron sheath between an emitter and a cesium plasma at equilibrium. 
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Figure 13. - Plane electron sheath between an emitter and a cesium plasma near equilibrium. 


59 



I II nil II I III I III III! Ill I 



X*E 4 

0. 

2S0. 

500. 

750. 

1000. 

1250. 

ISCC. 

1750. 

2000. 

2250. 

2500. 

0.C25 

O.C54 

C.116 

O.lSf 

0.195 

0.239 

0.287 

0.341 

C.40C 

0.466 

C- 1- 

1 
1 
1 
1 
1 
1 
1 
1 
1 

1 

1 

1 

1 

1 

1 

1 

1 

1 ♦ 

1 

1 

i 

1 

1 ♦ 

1 • ♦ 
*l ♦ 

* +l 

«• 1 

1 

1 

1 

* 

♦ 

1 

1 

1 

1 

♦ 1 ♦ 
* * 

* 1 

1 
1 
1 
1 
1 

1 • 

♦ 1 ♦ 
1 
1 
1 
1 
1 
1 
1 

* 4 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

*17 

0.544 

C.OC00501 

1 

1 ♦♦ 

1 

I 

1 

I 

1 

1 

I 

1 

i 

0.619 

1 

1 ♦ 

1 

1 

1 

1 

1 

1 

1 

1 


0.631 

1 

1* 

1 

1 

1 

1 

1 

1 

L 



0.721 

1 

* 

1 

1 

1 

1 

1 

1 

1 

1 

1 


0.134 

1 

*1 

1 

1 

1 

1 

1 





1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


0.856 

1 

♦♦1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

u 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1.C07 

C.CCClOOl * 

1 

1 

1 

1 

I 

I 

i 

1 

I 

1 

I 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

E 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

® 1.203 

1 * 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

E 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

o 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

<D 

o 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

™ 1-‘'72 

1 + 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

■K 1*^85 

C.CC01501 ♦ — 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


i 

1 

i 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

I 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1.S9C 

C .CCC2001 ♦ — 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

Effective 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

sheath ^^ 3 ^, 

1 

1 * 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

width = 2^32 7 

1 * 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 


1 

1 

1 

1 


1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

c. 

250. 

500. 

750. 

1000. 

125C. 

15C0. 

175C. 

2000. 

2250. 

2500X1 CAI 

NE1X10**-11) vs X 

* RICHAROSON-OUSHHXN 

* SCHOTTKY 




Electron number density, cm’^ 







Effective 
sheath 
width = 


-lUOOO. -9000. -8000. 

O.C25 C. 1 1 1 

C.C54 111 

O.llB 111 

0.15S 111 

0.19S 111 

0.2iS 111 

0.287 111 

0.341 111 

0.40C 1 1 1 

0.468 1 1 i 

C. 0000 501 1 1 

C.544 111 

0.619 111 

0.631 111 

C.721 111 

0.734 1 1 1 

1 1 I 

0.656 1 1 1 

1 1 1 

1 1 1 

1.007 C.OCOlOOl 1 1 

1 1 1 

1 1 1 

1 1 1 

1.203 1 1 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

1.472 111 

1.485 C.0CC1501 1 1 

1 1 1 

1 I 1 

1 1 1 

1 1 1 

1 1 1 

1 1 1 

111 
I 1 1 

1 1 1 

2.S9C C.CCC20C1 1 1 — 

1 1 1 

111 
1 1 1 

1 1 1 

1 1 1 

2.317 1+ 1 1 

2.327 1* 1 1 

1 1 1 

0.0002 50X10^1 1 1 

-10000. -9C00. -8000. 


hUXlO** -91 VS X 
♦ RICHAPOSON-OUShMXS 


♦ SCHOTTKY 


-7000. 


-6000. -5000. -40CC. -3000. -2000. -1000. 0. 

1 1 1 I 1 1 1 — 

1 1 1 1 I I *+ 1 

11111 1***+ 1 
1 1 1 1 1 I 

1 1 1 1 1 *♦! 1 

1 1 1 I 1 ♦ ♦ 1 1 

1 1 1 1 1 ♦ ♦ 1 1 

1 1 1 1 ♦ ♦ 1 1 

1 1 1 1 ♦ 1 1 
1 1 1 1 ♦ ♦ 1 1 1 

1 1 1 1 1 1 
1 1 1*1 1 1 1 

1 1 1*1 1 1 1 

1 1 14 1 1 1 1 

1 1*1 I 1 1 1 

1 1 1 1 1 I 1 

1 1 ♦ ♦ 1 1 1 1 1 

1 1 1 1 I 1 1 

1 1 1 I 1 1 1 

1 *— ♦— 1 1 1 1 1 1 — 

1111111 
1 1 1 1 1 I 1 

1 1 1 1 1 1 1 

* + 1 1 1 1 I 1 1_ 

1 1 1 I I 1 T 

1111111 
1111111 
1111111 
1 1 1 1 I 1 1 

i 1 1 1 1 1 1 

1 1 1 1 1 1 I 

1 1 1 1 1 I 1 

I 1 1 I I I 1 

1111111 
1111111 
1 I 1 1 1 1 1 

1111111 
1 1 1 1 1 1 I 

1 1 1 i 1 1 1 

1 1 1 1 1 1 I 

1111111 
1111111 
1 1 I 1 1 1 1 

1 1 1 1 1 I 1 

1 1 1 1 I 1 1 

1 1 1 I I I 1 

1 1 1 1 1 I 1 

-6000. -5CCC. -4CCC. -3CCC. -2000. -1000. OXIO^ 

Ion number density, cm‘3 


Figure 13. - Continued. 
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